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Glucagon, a peptide hormone secreted from the a-cells of 
the pancreatic islets, is critical for blood glucose homeostasis. 
We reviewed the literature and employed a computational 
systems analysis of intracellular metabolic and electrical 
regulation of glucagon secretion to better understand these 
processes. The mathematical model of a-cell metabolic 
parameters is based on our previous model for pancreatic 
|3-cells. We also formulated an ionic model for action 
potentials that incorporates Ca 2+ , K + , Na + and CI" currents. 
Metabolic and ionic models are coupled to the equations 
describing Ca 2+ homeostasis and glucagon secretion that 
depends on activation of specific voltage-gated Ca 2+ channels. 
Paracrine and endocrine regulations were analyzed with an 
emphasis on their effects on a hyperpolarization of membrane 
potential. This general model simulates and gives insight into 
the mechanisms of regulation of glucagon secretion under a 
wide range of experimental conditions. We also reviewed and 
analyzed dysfunctional mechanisms in a-cells to determine key 
pharmacological targets for modulating glucagon secretion in 
type 1 and 2 diabetes. 



1. Introduction 

Glucagon is the principal blood-glucose-increasing (counter- 
regulatory) hormone of mammals, with a variety of other 
effects as well. It is released from a-cells of the pancreatic islets, 
comprising only a few percent of the total cell mass of an islet. 
Electrophysiological studies on individual a-cells have revealed 
that they are electrically excitable. Stimulation of glucagon secre- 
tion is secondary to increased action potential (AP) firing and 
also may be affected by paracrine signals (mediated by factors 
released from neighboring (3- and 8-cells) including insulin, 
Zn 2t , gamma amino butyric acid (GABA) and somatostatin, 
by endocrine signals such as glucagon-like peptide-1 (GLP-1) 
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or leptin and by autonomic input (for a review, see refs. 1-5). 
Glucagon plays a critical role in glucose homeostasis, diabetic 
ketoacidosis and its secretion is dysregulated in hypoglycemia 
associated with diabetes. 4,6 Thus, it is important to elucidate the 
signals that trigger glucagon secretion and the transduction of 
these signals within the a-cell. 

Pancreatic a- and (3-cells show a completely different behav- 
ior depending on the serum glucose level. However, an inhibi- 
tory effect of glucose on the a-cell, like its stimulatory effect on 
the (3-cell, requires metabolism of the sugar. Indeed, the a-cells 
are equipped with many of the proteins involved in the meta- 
bolic glucose sensing of the (3-cell, including glucokinase 7 and 
ATP-regulated K + -channels (K ATp -channels). 2,8 Glucose may also 
increase the intracellular concentration of ATP in a-cells.'' 10 

Electrical activity in a-cells is difficult to study within an 
intact islet and so this process is typically studied as isolated puri- 
fied cells in vitro. Glucose concentrations modulate pancreatic 
a-cell electrical activity, Ca 2+ signal and hormone secretion. At 
low glucose concentrations (3 mM or below) isolated a-cells can 
be electrically active and AP firing is required for glucagon secre- 
tion. Increased glucose concentration paradoxically suppresses 
electrical activity whereas activity in (3-cells in increased. 1,11 ' 12 
This electrical activity usually involves a complex interaction 
between several different ion channels as in pancreatic (3-cells. 
Pancreatic a-cells express voltage-gated Na*, K* and Ca 2+ chan- 
nels that regulate AP firing at low glucose level. 1 ' 2 ' 8 

Despite the recent increase in our knowledge of a-cell physiol- 
ogy and biochemistry, we still lack a coherent model of regula- 
tion in a-cells. Significant discrepancies exist among studies of 
glucagon secretion. It is at present unclear how increased glucose 
metabolism is coupled to the suppression of electrical activity in 
the a-cell, and the topic is subject to conflicting reports with 
the possibility of species differences (for review, see refs. 1 and 
2). The present work attempts to bridge the gap between elec- 
trophysiological recordings, metabolic measurements, Ca 2+ han- 
dling and glucagon secretion. 

To achieve a better understanding of a-cell physiology and 
its role in diabetes we undertook an integrated consideration of 
data obtained on isolated a-cells, human and rodent islets and in 
clinical investigations. We used a computational systems biology 
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approach with construction of a mathematical model of intra- 
cellular processes. We have previously focused on applying this 
approach to mechanisms of glucose sensing of insulin secretion, 
Ca 2t dynamics, electrophysiology events and exocytosis in pan- 
creatic (3-cells. 13 " 18 Here we focus on pancreatic a-cells based on 
our computational models for (3-cells. This allows also a com- 
parison of the mechanisms of glucose sensing in a— and (3-cells. 
Our aim was the investigation of the properties of rodent and 
human a— cells with a help of a novel computational model for 
metabolic processes, electrical activity, Ca 2t handling and gluca- 
gon secretion. 

The mathematical model for a-cells was tested in various 
scenarios and the results of simulations were compared with the 
experimental data. The model gives a reasonable fit to important 
aspects of experimentally measured metabolic and AP proper- 
ties of a— cells and provides a framework for analyzing the role 
of metabolic changes and individual ionic currents on glucagon 
secretion. In additions to a direct effect of glucose, the influence 
of several paracrine and endocrine mechanisms on AP firing and 
glucagon release were considered. 

Defective regulation of glucagon secretion clearly contributes 
to type 1 and 2 diabetes in a variety of ways. 6,19 ' 20 For this reason 
there is significant interest in understanding how a-cells regu- 
late glucagon secretion in humans. We have attempted to analyze 
these processes and evaluate how particular impairments in the 
mechanisms of a-cell glucose sensing can lead to glucagon release 
changes in diabetes. We also discuss examples of potential thera- 
peutic strategies. 

2. Components of Glucagon Secretion Regulatory 
Mechanism and their Mathematical Description 

Here we will briefly outline the different parts of the glucagon 
secretion regulatory mechanisms and describe the equations and 
parameters utilizing in the simulations of a-cell function. A sche- 
matic diagram of the biochemical steps, Ca 2+ handling, chan- 
nels and mechanisms of glucagon secretion regulation in a-cell 
is presented in Figure 1. The data used to fit the computational 
model in our study were taken primarily from isolated rodent 
a-cells and islets. Since pancreatic a-cells remain incompletely 
understood, there remain gaps in any completely a-cell specific 
data set. For the purpose of making a convenient model we have 
chosen plausible suggestions based on mechanisms of regulation 
and parameters from (3-cells or other cell types to close the gaps 
pending additional data. 

2.1. Intrinsic metabolic regulation. The cellular metabolic 
mechanisms leading to insulin secretion in pancreatic (3-cells are 
fairly well understood (for reviews, see refs. 2 and 15). Similar 
machinery may operate in a-cells. Although both cell types pos- 
sess different glucose transporters, it has been demonstrated that 
glucose transport is not a limiting factor in the metabolism of 
this sugar. 21,22 Glucokinase was found in both cell types 7 and may 
be also the rate-limiting step for glycolysis in a-cells. 23 a-cells 
are equipped with ATP-sensitive K + channels of the same type as 
those constituting the resting conductance in (3-cells (see below 
in Sec. 3.3.2). It seems that a metabolic mechanism of sensitivity 



to glucose is similar in a- and (3-cells: mitochondrial metabolism 
generates ATP leading to an increase of the ATP/ADP concen- 
tration ratio, leading to closure of K ATp channels. The resulting 
decrease in K* efflux causes PM depolarization which changes the 
rate of exocytosis of glucagon granules. 2,24 

We have recently developed a mathematical model for reg- 
ulated glucose sensing in pancreatic (3-cells 15,25 and use it here 
to analyze a-cells. The cytoplasmic part of the model includes 
equations describing glucokinase, glycolysis, lactate dehydroge- 
nase, NAD(P)H and ATP production and consumption (Fig. 1). 
The mitochondrial part includes production of NADH, which is 
regulated by pyruvate dehydrogenase. NADH serves as a source 
of electron transfer in the chain of oxidative reactions to establish 
a proton motive force, driving the F ] F () ATPase. Redox shuttles 
and mitochondrial Ca 2 * handling were also modeled. We have 
performed a computational analysis of (3-cell fuel metabolism 
to quantitatively assess how cytoplasmic ATP/ADP ratio can be 
controlled by mitochondrial and cytoplasmic processes. 15,25 

Key differences in the rates of energetic processes were found 
in a-cells vs. (3-cells. Increased ATP/ADP ratios and NAD(P)H 
signals were found in rodent a-cells compared with (3-cells at low 
glucose levels. 23,26 In our previous (3-cell model of glucose sensitiv- 
ity no ATP production was generated without glucose consump- 
tion because only glucose served as a source of ATP production 
during glycolysis and NADH for oxidative phosphorylation. To 
describe a nonglycolitic source for ATP production in the absence 
of glucose or in low glucose we introduced a special term (other 
than pyruvate influx) in the equation for mitochondrial NADH 
([NADH] m ) dynamics in our model of glucose sensing. Now, the 
following equation describes the change in [NADH] m over time: 



d [NADH] /dt 
IAnHm [NADH] 



(4.6/„ YR +/ 



-0.1 L 



+ k )/V . - 



(Eqn. 1) 



where k ng is the novel coefficient describing a nonglycolitic source 
of NADH in mitochondria, / pYR is the rate of pyruvate consump- 
tion in mitochondria, / TNADH is the flux through the NADH 
shuttles, / hres is the rate of proton pumping through ETC, 4.6 
and 0.1 are the stoichiometric coefficients, V . is the relative 

mmir 

mitochondrial matrix fraction of the cell volume, k...„ u is the 
rate constant of NADH consumption in mitochondria (see for 
detail Eqn. 14 in ref. 15). 

Introducing Equation 1 instead of the corresponding equa- 
tion for NADP dynamics (Eqn. 14 in the model 15 ) we obtained a 
"complex" model for glucose sensing that can be used for a-cells 
and (3-cells simultaneously even at low glucose levels (see below). 
This model is distinguished from the previous model 15 only by k ng 
coefficient in Equation 1 and at k =0 these models are identical. 

2.2. Channels and plasma membrane potential regulation. 
Glucagon release, like insulin release, is influenced by physio- 
logical a-cell electrical activity that fundamentally resembles the 
excitation-secretion coupling seen in many secretory cell types 
(for reviews, see refs. 3 and 27). Like insulin-secreting (3-cells, 
a-cells express different types of ion channels. There are at least 
four different types of IC-selective channels: the ATP-sensitive 
K + -channel (K ATp -channel), the delayed rectifying K + -channel 
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(K Dr -channel), a transient K + -channel 
(A-type K + channel) and a G protein-gated 
K + -channel activated by somatostatin 
(K^-channel). Other channels include 
voltage-gated Na + -channels, at least four 
types of voltage-gated Ca 2+ -channel 
(VGCC) (T-, N-, P/Q- and L-type) and 
GABA A receptor chloride-channels (refs. 
1,2,8,28; see also below Sec. 3.3. for 
details). In the Supplemental Material we 
describe the properties (mostly similar to 
that described in other cells) of each type 
of ion channel and their mathematical 
descriptions. 

A schematic diagram of the a-cell spe- 
cific channels is presented (Fig. 1). The 
AP in a single a-cell can be described 
with the following current balance differ- 
ential equation 



- C m ( dV P /dt ) = 7 KATP + 7 KD, + 7 KA + 
7 KGS + 7 CaL + 7 CaNL + 7 CaT + 7 PCa + 7 Na 



ck; 



+ / Nab 

(Eqn. 2) 



where V is the plasma membrane 
potential, t is time, C m is the whole-cell 
membrane capacitance, / KATP is the ATP- 
sensitive K* channel current, I KDi is the 
rapid delayed-rectifier K + current, 7 KA 
is the high voltage-activated A-type K + 
current, I KGS is the K* current through 
G protein-gated K*-channel activated 
by somatostatin, 7 CaL is the high-voltage- 
activated L-type Ca 2t current, / CaNL is the 
high-voltage-activated non-L-type Ca 2+ 
currents, / CaT is the low-voltage-activated 
"T-type" Ca 2t current, / pCa is the plasma 
membrane Ca 2t -pump current, 7 Na is the 
voltage-gated Na* current, 7 Nab is the Na + 
background current, 7 C|G is the GABA 
activated CI' current. 

2.3. Ca 2+ handling. Ca 2t enters 
into cells through Ca 2+ channels and is 
removed by Ca 2+ pumps and exchangers. 
Ca 2t handling can also include ER and 
mitochondrial sequestration. Including 
only fluxes through Ca 2t channels and 
Ca 2t pumps on the PM, the equations for 
free cytoplasmic Ca 2+ ([Ca 2+ ] c ) dynamics 
can be written: 



d[Ca 2 1 c /d.= f (-/ CaL -/ CaNL -/ CaT -2 

(Eqn. 3) 
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Figure 1. Schematic diagram of the ionic current, metabolic processes, Ca 2+ fluxes and glucagon 
secretion mechanisms in pancreatic a-cells. Transmembrane currents are: / KATp is the ATP-sensitive 
K + channel current, / KDr is the rapid delayed-rectifier K + current, / is the high voltage-activated 
A-type K + current, / KGS is the K + current through G protein K + -channel activated by somatostatin, 
/ CaL and / CaNL and are high-voltage-activated L-type and non L-type Ca 2+ currents, / CaT is the low- 
voltage-activated "T-type" Ca 2+ current, / pCa is the plasma membrane Ca 2+ -pump current, / is the 
voltage-gated Na + current, / Nab is the Na + background current, / is the GABA activated CI cur- 
rent. Glucose equilibrates across the plasma membrane and is phosphorylated by glucokinase to 
glucose 6-phosphate, which initiates glycolysis. Lactate dehydrogenase (LDH) converts a portion 
of pyruvate to lactate. Pyruvate produced by glycolysis enters the mitochondria and is metabo- 
lized in the tricarboxylic acid (TCA) cycle, which then yields reducing equivalents in the form of 
NADH and FADH2. k ng is the coefficient describing a nonglycolitic source of NADH in mitochondria. 
The transfer of electrons from these reducing equivalents through the mitochondrial electron 
transport chain is coupled with the pumping of protons from the mitochondrial matrix to the 
intermembrane space. The resulting transmembrane electrochemical gradient drives the ATP 
synthesis at ATP-synthase. Part of the protons may leak back through uncoupling proteins (UCPs). 
The shuttle systems are required for the transfer of reducing equivalents from the cytoplasm to 
the mitochondrial matrix. ATP is transferred to the cytosol, raising the ATP/ADP ratio. This results 
in the closure of the ATP sensitive K + channels, which in turn leads to depolarization of the cell 
membrane. In response, the voltage-sensitive Ca 2+ channels open, promoting calcium entry and 
triggering exocytosis of glucagon granules. ATP c and ADP free are the free cytosolic form of ATP 
and ADP, G3P is the glyceraldehydes 3-phosphate, PDH is the pyruvate dehydrogenase, ANT is the 
adenine nucleotide translocase, ¥ is the mitochondrial membrane potential. Solid lines indicate 
flux of substrates, and dashed lines indicate regulating effects, where (+) represents activation and 
(-) repression. 



where f is the fraction of free Ca 2+ in 
cytoplasm, F is Faraday's constant, V. is 
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Table 1. Cell parameters 



Parameter 


Definition 


Value 


Equation 


Reference 


r 


Membrane capacitance 


5,000 fF 


2 


a 


v 


Cytosol volume 


0.764 pL 


3 


b 


f 


Fraction of free Ca 2+ in cytoplasm 


0.01 Ul 


3 


b 


F 


Faraday constant 


96,487.0 C/mol 


3 




k 

sg 


Coefficient of the sequestration rate of [Ca 2+ ] 


0.0001 ms _1 


3 


b 


N 

a 


Amount of a-cells in one islet 


200 


4 


f 


V 

pg 


Volume of the vessel 


1000 Ul 


4 


f 


k 

eg 
k 

rg 


Stoichiometric coefficient 


1.4 Ul 


5 


f 


Glucagon production in rest conditions 


0.0003 Ul 


5 


f 




Perfusion coefficient 


0.0005 ms 1 


6 


f 


k 

pg 


Scaling coefficient 


2000 Ul 


6 


f 



a 44 ; b 17 ; f, adjusted to fit the experimental values; Ul, unitless. 

the effective volumes of the cytosolic compartment, and k is a 
coefficient of the sequestration rate of Ca 2t at exocytosis. 

However, detailed data on Ca 2t dynamics during one spike 
in a-cells is lacking, so we used the coefficients we employed for 
mouse (3-cells electrophysiological model 14 (Table 1). 

2.4. Role of Ca 2+ channels in glucagon exocytosis. 
Capacitance measurements in a-cells have shown that the glu- 
cagon granule exocytosis is negligible at membrane potentials 
more negative than -20 mV but increases steeply at more depo- 
larized voltages. 3 ' 8 ' 29 ' 30 This behavior may reflect a necessity for 
local calcium entry through high voltage-gated non-L-type Ca 2+ 
channels in rodent and human a-cells because the process of 
release of primed granules in a-cells (as well as in (3-cells) may 
need a local highly increased calcium in specific microdomains 
(for a review, see refs. 29 and 31). Microdomain Ca 2+ may be 
proportional to the Ca 2+ current through local Ca 2t channel, 32 
and we model granule exocytosis to a first approximation as the 
dependence on Ca 2+ influx through specific non-L-type Ca 2+ 
channels. 

Here we modeled the electrophysiological events inside 
one cell. However, in a single (3-ceIl a few insulin granules are 
released per minute 16,33 and a-cells are likely to have similar char- 
acteristics. 5,34 At present it is difficult to model a-cell exocytosis 
in single cells (see below, Sec. 4.5). Measured rates of glucagon 
secretion are also very variable and often represent only a com- 
parison of relative rates. However, we found that an evaluation of 
the relative rates of glucagon secretion can be sufficient for our 
consideration. 

For these reasons, we have adopted the convention of simu- 
lating one islet containing many a-cells with which to simulate 
a relative rate of glucagon secretion. We simulated a hypotheti- 
cal perfusion experiment, where one islet is localized in a perfu- 
sion chamber. In this case the relative glucagon concentration in 
media surrounding this islet and relative glucagon output from 
the chamber can be written as 



(-/„ NI /F) + * 



G = k k [Glc] 



(Eqn. 5) 
(Eqn. 6) 



d[Glc] / dt = f N / V -k [Glc] 

J s gs a pg pc J s 

where 



(Eqn. 4) 



[Glc] is the relative glucagon concentration in medium, N a is 
the amount of a-cells in one islet, f is the rate of exocytosis in 
one cell, V is the volume of the chamber, k is the stoichiometric 
coefficient, / CaNL /F is the Ca 2+ influx through specific VGCCs 
(non-L-type Ca 2+ channels in this work), is the rate coefficient 
for relative glucagon production in rest conditions (in hyperpo- 
larized cells), G s is the relative glucagon secretion rate in the per- 
fusion experiment, k is the perfusion coefficient and k is the 
r ' p C r p g 

scaling coefficient. 

To compare with the experimentally measured rates of glu- 
cagon secretion we also propose that a-cells in an islet account 
for approximately 20% of islet cells, i.e., about 200 a-cells in a 
mouse islet if an islet has about 1000 cells 33 (coefficient N a in 
Eqn. 4). Other unknown coefficients in Equations 4-6 were 
evaluated according to the proposal that the mean relative gluca- 
gon concentration in the medium and the relative secretion rate 
were about 1 relative unit at low glucose in our hypothetical per- 
fusion experiment (Table 1). 

2.5. Computational model of regulation of glucagon secre- 
tion. In parallel with the metabolic model of glucose sensing, we 
have constructed a computational model of glucagon secretion in 
a-cells as the combination of the equations for the regulation of 
PM potential (Eqn. 2), Ca 2+ handling (Eqn. 3), relative glucagon 
release (Eqns. 4-6) and the voltage-dependent gating variables 
for channels (see Supp. Material). The influence of glucose was 
modeled as the changes in ATP/ADP ratio obtained from the 
model of glucose sensing that determines K AT[) channel opening. 

However, precisely detailed mathematical models of electro- 
physiological, metabolic and exocytotic processes are impossible 
to construct at the present time (and in the immediate future) due 
to lack of both adequate experimental data and detailed models 
of separate processes. We also evaluated the model parameters 
from the literature, when possible, but they were also found by 
fitting a set of known experimental data as an important reality 
test (see Sees. 2.1-2.4 and Supplemental Material). 
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Table 2. Membrane current parameters 
Parameter 



9mKATP 
9 m KDr 
T dKr 
^dKr 
k dKr 
9 m KA 

v d , 
K, 

T fK 
^fKA 
^fKA 
9mKGS 

E„. 



dKA 
dKA 
KiKA 
T KA 



"f2Na 

9mN, 
E, 



ab 



=>mCaL 



"fCaL 
KcaL 
9mCaNL 
T dCaNL 
^dNL 



Definition 

Reversal potential for K + current 
Maximum conductance for / KATp 
Maximum conductance for / DKr 
Time constant for d Kr 
Half-activation potential 
Slope of half-activation potential 
Maximum conductance for / 
Time constant for d Kfl 
Half-inactivation potential 
Slope of half-activation potential 
Time constant for f 
Half-inactivation potential 
Slope of half-inactivation potential 
Maximum conductance for / KGS 

Reversal potential for Na + 
Maximum conductance for /., 

Na 

Time constant for d„ 

Na 

Half-activation potential 
Slope of half-activation potential 
Time constant for f,.. 

INa 

Half-inactivation potential 
Slope of half-inactivation potential 
Time constant for f,„ 

2Na 

Half-inactivation potential 
Slope of half-in activation potential 

Maximum conductance for / Nab 
Reversal potential for Ca 2+ current 
Conductance for L , 

CaL 

Half-activation potential 
Slope of half-activation potential 
Time constant for f,. , 

ICaL 

Half-inactivation potential 
Slope of half-inactivation potential 
Conductance for / CaNL 
Time constant for d CaNL 
Half-activation potential 
Slope of half-activation potential 
Time constant for f CaNL 
Half-inactivation potential 
Slope of half-inactivation potential 
Maximum conductance for /. . 

CaT 

Time constant for d r _ 

CaT 

Half-activation potential 
Slope of half-activation potential 
Time constant for f r _ 

CaT 

Half-inactivation potential 



Value 

-75 mV 
30,000 pS 
18,000 pS 
18 ms 
-15 mV 
7 mV 
1200 pS 
0.26 ms 
-9 mV 
8mV 
12 ms 
-20 mV 
4.8 mV 
OpS 
70 mV 
32000 
0.15 ms 
-30 mV 
10 mV 
1 ms 
-42 mV 
6mV 
40 ms 
-44 mV 
6.0 mV 
44 pS 
70 mV 
100 pS 
-20 mV 
6mV 
6.8 ms 
-20 mV 
6.5 mV 
300 pS 
0.4 ms 
-9 mV 
8mV 
18 ms 
-9 mV 
8mV 
200 pS 
0.4 ms 
-50 mV 
6mV 
6.8 ms 
-74 mV 



Equation 

A1 

A4 
A5 
A6 

A8 
A9 
A10 

A11 
A12 

A13 
A14 
A14 
A15 
A16 

A17 
A18 

A19 
A20 

A21 

A22 

A24 

A26 
A27 

A28 
A29 
A30 

A31 
A32 

A33 
A34 
A35 

A36 
A37 



Reference 

a 



a"; b 8 ; f, adjusted to fit the experimental values. 
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Table 2. Membrane current parameters (continued) 



Parameter 


Definition 


Value 


Equation 


Reference 


k fCaT 


Slope of half-inactivation potential 


6mV 




f 


P . 
mCap 


Maximum / pCa current 


300 fA 


A38 


f 


^Cap 


Half-maximum Ca 2+ binding constant for / pCa 


0.3 uM 




f 


9mCIG 


Maximum conductance for / 


OpS 


A39 


f 


Ec, 


Reversal potential for CI" 


-60 mV 


A39 


f 



a"; b 8 ; f, adjusted to fit the experimental values. 
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Figure 2. Effect of increasing glucose on a- and (3-cells energetics. Ex- 
tracellular glucose concentration was varied and the steady-state simu- 
lations of the model parameters are represented. Simulations were run 
with the basic set of parameters for (3-cells (Tables 2 and 3 from ref. 15) 
and additionally with the coefficient describing the nonglycolitic source 
of NADH in mitochondria (k ng , Eqn. 1). Corresponding parameters were 
changed for a-cells: k n9 was increased from 0.006 (jimol ms 1 (for (3-cells) 
to 0.05 txmol ms V mLDH (maximal lactate dehydrogenase activity) was 
increased from 1.2 (jimol ms" 1 to 12 (jimol ms 1 , the maximal pyruvate 
dehydrogenase (PDH) activity (V mpdh ) was decreased from 0.3 (jimol ms ' 
to 0.08 (jimol ms 1 , the regulated leak coefficient (P |r ) was increased from 
0.0012 (jimol ms 1 to 0.008 (jimol ms 1 , the glucokinase activity (V J was 

taken as 0.01 (aitioI ms 1 for both a- and (3-cells. ( ) represents (3- and 

( ) a-cells. (A) (••••) is the rate of the glucokinase reaction (close in 

a- and (3-cells), J pYR is the rate of pyruvate consumption in mitochondria 
for (3- and a-cells. (B) M" 1 is the mitochondrial membrane potential. (C) 
[NADH] m is mitochondrial NADH. (D) [ATP] c /[ADP] c is the ATP/ADP ratio 
in the cytoplasm. 



Tables 1 and 2 contain all the information necessary to carry 
out the simulations presented in this paper and are pointed out 
in the text as a simulation at basal level. In some cases, when the 
values of parameters differed from those at basal levels the values 
assigned are noted in the text or in the corresponding figure leg- 
ends. To calculate the steady-state cellular parameters, the model 
was allowed to run up to steady values with no changes in coef- 
ficients. The units, except where indicated otherwise, are: time in 
milliseconds (ms), voltage in millivolts (mV), concentration in 
micromoles /liter (|xM), current in femtoamperes (fA), fluxes in 
mmol s" 1 , conductance in picosiemens (pS), capacitance in femto- 
farads (fiF). Numerical integration was performed using standard 
numerical methods. 

These models are available for direct simulation on the web- 
site "Virtual Cell" (www.nrcam.uchc.edu) in "MathModel 
Database" on the "math workspace" in the library "Fridlyand" 
with the name: "Alpha and beta cells glucose sensing" and 
"Glucagon secretion" are represented. Visualization and graphi- 
cal analysis were performed using "Excel." 

3. Review, Simulations and Analysis of Experimental 
Observations of a-cell Intracellular Regulatory 
Mechanisms 

In this section we use computational simulations to critically 
review and analyze the experimental observations of intracellu- 
lar regulatory mechanisms, where electrophysiological data were 
obtained for isolated a-cells. However, because our mathematical 
model of glucagon secretion has significant limitations (see Sec. 
2.5) we will use it mainly to test and explain numerous seemingly 
mutually incompatible experimental observations and the devel- 
opment of hypotheses concerning the mechanisms of regulation of 
glucagon secretion, not for simulation of individual experiments. 

Here, we define a "low glucose concentration" as approxi- 
mately 1—2 mM (where glucagon secretion is maximal in isolated 
a-cells), "medium glucose concentration" is 3-5 mM (where glu- 
cagon secretion is minimal) and "high glucose concentration" is 
above 10 mM (where glucagon secretion begins to increase) (see 
below in Sec. 3.2). 

3.1. Intrinsic metabolic regulation in a-cells. The ratio of 
lactate dehydrogenase to mitochondrial glycerol phosphate dehy- 
drogenase activity is higher in a- compared with (3-cells. 35,36 (Low 
activity of lactate dehydrogenase is favoring the mitochondrial 
oxidation of glucose in (3-cells 15,35 ). In contrast, pyruvate carbox- 
ylase, that is required for conversion of pyruvate to oxaloacetate 
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for repleting the Krebs cycle with intermediates, has significantly 
higher expression in (3-cells. 36 Interestingly in this regard, uncou- 
pling protein 2 is expressed much more highly in a-cells com- 
pared with (3-cells. 37 

We examined how these differences can affect the meta- 
bolic responses of a-and (3-cells to glucose, using the corrected 
model for glucose sensing (Fig. 2). To describe characteristic 
peculiarities of a-cells the activity of pyruvate dehydrogenase 
was decreased, and the activities of lactate dehydrogenase and 
uncoupling protein 2 were increased to compare with the (3-cell 
model. 15 The rate of NADH production from nonglucose sources 
was assumed to be higher in a- compared with (3-cells to take 
into account the increased ATP and NAD(P)H levels in a-cells 
at low glucose levels (see Sec. 2.1). We have found that in this 
case the rate of pyruvate consumption in mitochondria (/ PYR ) was 
significantly decreased in a-cells (Fig. 2A). Mitochondrial mem- 
brane potential, ATP/ADP ratio and NADH in mitochondria 
were higher at low glucose and saturated sharply in ct-cells unlike 
(3-cells simulation (Fig. 2B-D). 

The results of simulations are consistent with experimental 
data: glycolytic flux, measured as D-[5 3 H] glucose utilization, 
was comparable in a- and (3-cells. Rates of glucose utilization 
increased proportionally to the extracellular substrate concentra- 
tion between 1 and 10 mM substrate and were virtually iden- 
tical in non-(3-cells (mainly a-cells) and (3-cells. 22 The ATP/ 
ADP ratio was significantly higher in a-cells compared with 
(3-cells at low glucose levels. 26 For both a- and (3-cells, increased 
glucose concentration leads to an elevation of the intracellular 
ATP/ADP-ratio. 1,2,38 However, experimental data on glucose 
metabolism agrees with the observation that ATP/ADP changes 
in response to glucose are less variable in a-cells compared with 
(3-cells. 9,10,26 The NAD(P)H signal was greater in mouse islet 
a-cells compared with (3-cells at 1 mM glucose, however, the 
increase of this signal with glucose was less in a-cells. 23 Purified 
rat non-(3 cells (approximately 80% glucagon-producing a-cells) 
exhibit a much lower rate of glucose oxidation (only 30%) in 
response to glucose than in purified (3-cells, 36 corresponding to 
the decreased rate of pyruvate consumption in mitochondria in 
our model (Fig. 2A). Therefore, our systems analysis indicates 
that recorded differences in metabolic response between a-cells 
and (3-cells can be explained in part by the variations in expres- 
sion of the specific enzymes that are responsible for the glucose 
sensitivity mechanism. Based on this evaluation we can suggest 
that a-cells might have significantly suppressed mechanism of 
K AT[) channel inhibition with increased ATP/ADP ratio in com- 
parison with (3-cells. 

3.2. Spontaneous AP firing in a-cells and metabolic mecha- 
nisms of regulation of glucagon secretion. Patch-clamp experi- 
ments on isolated mouse, rat, guinea pig and human a-cells have 
show n that, unlike (3-cells, they produce spontaneous AP firing 
in the absence or low glucose levels (see refs. 1, 8, 39 and 40). 
We were able to simulate AP firing at low glucose levels (at the 
low [ATP] /[ADP] ratio) using the specified basal set of chan- 
nels and coefficients (Fig. 3). Our model shows that the PM can 
be depolarized in a-cells up to a threshold for AP firing even 
with low glucose levels. Simulated AP spikes were in the range 
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Figure 3. Modeling of spontaneous spikes, glucagon secretion and 
the changes of intracellular Ca 2+ concentration in response to a step 
increase in glucose concentration in a-cells. (A) PM voltage (V p ) and cor- 
responding detailed spike picture. (B) Relative glucagon secretion rate 
(Gs). (C) Cytoplasmic calcium ([Ca 2+ ] c ) transient. Glucose-induced K flTp 
conductivity changes were simulated by a step increase of the [ATP]/ 
[ADP] ratio at the arrow (1) from a low (ATP/ADP = 2, that corresponds 
to 1.3 mM glucose in Fig. 2D) to an intermediate glucose level (ATP/ 
ADP = 6, that corresponds to 3.6 mM glucose in Figure 2D) and at the 
arrow (2) from an intermediate to high glucose level (ATP/ADP = 19, that 
corresponds to 25 mM glucose in Fig. 2D). All other parameters were 
taken from the basic set of parameters (Tables 1 and 2). Changes in V p 
were simulated using Equation 1, [Ca 2+ ] using Equation 3, and relative 
glucagon secretion using Equation 6. 



between -60 and 10 mV. Note, the model responds with depolar- 
ization and generation of AP firing in the same voltage range and 
frequency as was observed experimentally in mouse and human 
a-cells (see for example, refs. 8, 11 and 41). 

Isolated mouse and rat pancreatic a-cells were electrically 
active in the absence of glucose but glucose inhibited electrical 
activity and suppressed glucagon secretion. 24,40 We have also sim- 
ulated the increase of [ATP] /[ADP] ratios following a stepwise 
increase of glucose from low to medium and high levels (arrows 
1 and 2 in Fig. 3B). Increased [ATP] /[ADP] ratio resulted in the 
closure of K ATp channels similarly as in the (3-cell mathematical 
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Figure 4. Detailed simulated curves of currents during a single spike 
from Figure 3A. (A) Action potential (V p ); (B) Principal currents for 
action potential (/ , / and LJ are represented for one characteristic 
spike, where / is main depolarizing and / and / KDr are main repolar- 
izing factors. (C) Ca 2+ currents: / CaL and / CaNL . 



model (compare refs. 15 and 17) leading to PM depolarization. 
However, quite the opposite of the case in (3-cells, PM depolariza- 
tion can lead to suppression of existing AP firing in a-cells. This 
leads to reduced opening of the voltage-gated calcium channels, 
decreasing Ca 2+ influx through specific VGCCs (non-L-type in 
mouse) and suppression of glucagon granule exocytosis with the 
increase of glucose from low to medium levels. The simulations 
correspond to the experimental observations discussed above. 

We will denote this mechanism of glucose action as a 
metabolic mechanism of regulation of glucagon secretion. 
According to this mechanism decreased glucagon secretion with 
increased glucose occurs as a consequence of increased ATP/ 
ADP and block of K ATp channels leading to PM depolarization 
and suppression of AP firing that exists at low glucose level in 
a-cells. 

In addition to the inhibitory effect of glucose seen at medium 
(3—5 mM) levels, glucose also exerts a stimulatory effect on glu- 
cagon release at higher concentration. For example, glucagon 
secretion at 20 mM was significantly higher than that observed at 
3-5 mM glucose [but it was lower than at low glucose (1-2 mM)] 



in isolated mouse and human islets 



Although the precise 



mechanism of this stimulatory effect is still obscure, our simula- 
tion shows that an additional suppression of K ATp channels with 
increased glucose up to "high levels" results in PM depolarization 
and increased glucagon secretion. This may be a result of increas- 
ing voltage-dependent activation of non-L-type Ca 2+ channels 
even without activation of AP firing (Fig. 3). 

3.3. Channels as the modulators of PM potential, Ca 2+ han- 
dling and glucagon secretion. We have evaluated here only the 



role of the most well-studied currents in a-cells. Taking into 
account this set of channels and currents we can reasonably 
describe the existing a-cell electrophysiological and Ca 2 * han- 
dling properties together with glucagon secretion data. 

3.3.1. Mechanism of AP spikes. A sequence of steps describing 
the mechanism of AP spikes in a-cells during stimulus-secretion 
coupling has been proposed. 2,24,44 According to this suggestion, 
K ATp channels are responsible for changes in the membrane 
potential. At low glucose concentrations, the negative membrane 
potential in these cells is situated within a range that allows 
activation of T-type voltage-dependent Ca 2+ channels. These 
channels drive the membrane potential to more positive val- 
ues, allowing the activation of Na* and VGCCs. Repolarization 
would take place following activation of A- and K dr -type K* chan- 
nel. An increased extracellular glucose level would produce a rise 
in the ATP/ADP ratio as a consequence of glucose intracellular 
metabolism, and block K ATp channels. This inhibition would 
allow the PM to depolarize to potentials that change the entire 
panoply of currents sustaining the AP, consequently suppressing 
electrical activity. 

Our simulation confirms the possibility for this sequence of 
steps and allows a detailed analysis to be performed for the role 
of specific channels. The voltage-dependent currents during one 
cycle of spontaneous activity in the ionic model have been ana- 
lyzed (Fig. 4). The upstroke phase is clearly initially generated by 
voltage-dependent Na* channels (/ Na ). When Na + enters the cell, 
the membrane is further depolarized and high-voltage-activated 
Ca 2 * channels open and activate transiently. The AP repolar- 
ization phase initially includes an activation of fast A-type K + 
channels. After that the rapid delayed rectifier K + channels also 
activate (7 KDr current) leading to increased repolarization current. 
The repolarization phase also includes the effects of the time- 
dependent inactivation of Na* and Ca 2+ channels. This restores 
the resting potential and transiently hyperpolarises the a-cell 
(after-hyperpolarization) . The after-hyperpolarization allows 
complete reactivation of the Na*-channels and Ca*-channels so 
that the next action potential is of full size. 

Voltage-dependent Na* and K* channels (the delayed-rectifier 
and A-type) together may define the spike-generation mecha- 
nism in a-cells. Other currents included in the model were rela- 
tively small during a spike generation and had an insignificant 
influence on AP shape. This simulated pattern corresponds to 
the experimental observation that the Na*-current amplitude in 
a-cells in freshly isolated mouse islets is 5-10-fold larger than 
combined Ca* currents. 31 However, small currents can play an 
important role in total interspike current and influence on spike 
frequency and height (see below). 

Glucagon secretion rate depends on activation of specific 
high voltage-gated calcium channels (non-L-type, see discussion 
above in Sec. 2.4). Simulated Ca 2 * current through these chan- 
nels during one spike is represented in Figure 4C. (However, 
the T-type Ca* current was relative small and is not shown). A 
decrease in amplitude or complete inhibition of AP firing will 
therefore result in small openings of these specific VDCCs and 
this in turn means that fewer glucagon-containing granules 
undergo exocytosis. 
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ATp -channels (Fig. 5). Note, in 
channels do not 



3.3.2. K ATp channels. Glucose decteases K ATp channel conduc- 
tance through an increased ATP/ADP ratio that was simulated 
in our model (Fig. 3). However, AP firing in isolated human and 
rodent a-cells and glucagon release in intact rodent and human 
islets also can be suppressed both when K ATp -channels are maxi- 
mally activated by a high concentration of diazoxide (an activator 
of the K ATp -channel) or when they are inhibited by a high con- 
centration of tolbutamide (a blocker of the K ATp -channel) at low 
glucose levels. 3 ' 8 - 24 ' 27 ' 44 - 45 

We simulated these experimental conditions by changing the 
maximal conductivity of K ATp -channels and have obtained the 
results similarly to the experiments described above. That is, AP 
firing and glucagon secretion was suppressed both under condi- 
tions of inhibited and activated K 
our model based on the experimental data, K ATp 
open fully at 0 mM glucose and close at high glucose because the 
ATP/ADP ratio has a sufficiently high level even at low glucose 
(as discussed above in Sec. 3.1) and the ATP/ADP (ratio) increase 
cannot fully block the K ATp -channels (see Supp. Material, Eqns. 
Al— A3). For this reason, the inhibitory effect of diazoxide on 
glucagon secretion can be explained by an additional PM hyper- 
polarization (as a consequence of increased K ATp -channel conduc- 
tance) up to a level that is lower than the threshold for AP firing. 
Abolishing AP firing and the associated decreased glucagon 
secretion in a simulation of additional K ATp -channel inhibition 
(that corresponds to the action of tolbutamide) can be explained 
similarly to the action of increased glucose as a consequence of 
PM depolarization (see above in Sec. 3.2). 

3.3.3. Delayed rectifying K*-channel. Inhibition of K Dr - 
channels results in suppression of glucagon secretion at low, but 
not high, glucose levels in mouse and human isolated islets. 8,46 
We simulated these experiments by decreasing K Dr channel 
conductance that corresponds to the effect of specific blockers 
of these channels (for example, stromatoxin) (Fig. 6). As sug- 
gested 8 ' 46 a decrease in maximal K Dr channels conductance sup- 
pressed AP firing and decreased glucagon secretion. This occurs 
because the voltage-gated K Dr channels contribute significantly to 
the AP downstroke in a-cells (see above in Sec. 3.3.1) and after 
their inhibition there is no reactivation of the voltage-gated Na + 
and Ca 2t -channels. This paradoxically suppressed AP firing even 
though this leads to PM depolarization. At higher glucose levels 
AP firing is inhibited and block of K Dr channel did not affect 
glucagon secretion (simulation not shown). 

3.3.4. G protein-gated K*-cbannel and somatostatin action. 
Somatostatin inhibits glucagon secretion in pancreatic a-cells. 39 ' 47 " 49 
Electrophysiological recordings showed that somatostatin acti- 
vated a low-conductance G protein-coupled K*-channel in single 
rat and mouse a-cells, causing hyperpolarization of the PM and 
inhibition of electrical activity. 39,47,48 We have included this cur- 
rent in our model (I KGS see Fig. 1 and Supp. Material, Eqn. A13). 
This current acts on the PM potential similarly to the K ATp chan- 
nel current, because it can be described by similar equation (com- 
pare Eqns. Al and A13 in Supp. Material), that is, its activation 
can lead to PM hyperpolarization. Similarly to activation of K ATp 
channels this can result in inhibition of AP firing and a decrease 
in glucagon release. An example of a simulation of somatostatin 




10 20 30 40 
time (sec) 



Figure 5. Simulated parameters at activation and block of K flTp chan- 
nels. (A) PM voltage (V p ) and (B) glucagon secretion rate transient in 
response to the changes in K ATp channel conductance at low glucose 
level (ATP/ADP = 2). Initially the maximal conductance of K ATp channel 
(g mKflTP , Eqn. Al) was increased from basal level up to 60 nS to simulate 
K flTp channels activation and PM repolarization. Then g mKflTP was de- 
creased to basal level (30 nS, Table 2) at arrow 1 that leads to AP firing 
and an increase in relative glucagon secretion. For simulation of K flTP 
channel block g „,„ was decreased from 30 nS to 10 nS at arrow 2 that 
leads to PM depolarization, AP firing suppression and decreased rela- 
tive glucagon release. Then g mKATP was further decreased from 10 nS to 2 
nS at arrow 3 leading to additional PM depolarization with an increased 
glucagon secretion as a consequence of supplementary activation of 
non-L type Ca 2+ channels with increased PM potential. Other coeffi- 
cients and initial conditions were as in Figure 3. 
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Figure 6. Simulation of block of rapid delayed rectified voltage-gated 
K + (K Dr ) channels. (A) PM voltage (V p ) and (B) glucagon secretion rate 
transient in response to a decrease of K channel conductance at low 
glucose level (ATP/ADP = 2). For simulation of block of K Dr channel the 
maximal conductance (g mKr , Eqn. A4) was decreased from 18 nS (basal 
level) to 0 nS at arrow leading to PM depolarization, AP firing sup- 
pression and a decrease in relative glucagon secretion. Low glucose 
induced AP firing was simulated initially as in Figure 3. 
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Figure 7. Simulation of activation of G protein-gated K + -channel acti- 
vated by somatostatin. (A) PM potential (V p ) and (B) relative glucagon 
secretion (Gs). For simulation of G protein-gated K + -channel activa- 
tion the maximal conductance (g Eqn. A13) was increased from 
0 nS (basal level) to 0.5 nS at arrow leading to PM repolarization, AP 
firing suppression and a decrease in glucagon secretion. Low glucose 
induced AP firing was simulated initially as in Figure 3. 




Figure 8. Simulated glucose induced spikes behavior and glucagon se- 
cretion at block of Na + channel. (A) Low glucose induced AP firing (Vp) 
with the initial simulation as in Figure 3. (B) relative glucagon secretion 
rate transients. For simulation Na + channel block the maximal conduc- 
tance gmNa (Eqn. A14) was decreased from 32 nS (basal level, Table 2) 
to 0 nS at arrow that leads to AP firing suppression and a decrease in 
glucagon secretion. 



action is shown in Figure 7. This is consistent with the inhibitory 
effect of somatostatin on glucagon secretion. 

3.3.5. Voltage-gated Na* channels. Tetrodotoxin, a specific 
blocker of certain voltage-gated Na* channels, strongly inhibited 
glucagon secretion in mouse and human isolated islets at low 
but not high glucose levels. 8 ' 24,30 ' 44 We utilized the mathemati- 
cal model presented here to estimate the impact of Na + channels 
on AP firing and glucagon secretion (Fig. 8). Blockade of Na + 
channels induced partial PM repolarization. However, the main 



effect of Na* channel block was a suppression AP firing and con- 
sequently glucagon secretion decreased. This occurs because the 
voltage-gated Na*-channels are critical for the development of the 
depolarizing phase of AP in a-cell (see above in Sec. 3.3.1). At 
high glucose levels AP firing was still inhibited and Na* channel 
block did not change glucagon secretion significantly (not shown) . 

3.3.6. Na* background current. Closure of K ATp -channels alone 
is not sufficient to cause membrane depolarization. Therefore, 
(3-cells must be equipped with some sort of depolarizing inward 
current. However, the nature of the depolarizing current remains 
elusive. 50 For example this current in (3-cells might result from for- 
ward operation of the Na-Ca exchanger, i.e., this can cause a Na* 
inward current because one Ca 2 * from the cytoplasm exchanges 
for three Na* from the medium." Another suggestion was activa- 
tion of cationic (principally Na*) 7r^-channels. 52 ' 53 In this article 
we used a model of some generalized inward Na* current (with- 
out specifying exactly which one) to simulate this depolarizing 
current in ot-cells (Eqn. 21 in Supp.Material). However, clearly 
additional experiments are required to define the properties of 
this current and its influence on a-cells in greater detail. 

3.3.7. L-type Ca 2 * channels. Glucagon secretion evoked by low 
glucose alone is "almost resistant" or only slightly affected by 
L-type Ca 2 * channel blockers in isolated mouse islets. 30 Decreased 
L-type Ca 2 * channel conductance also did not suppress AP firing 
and glucagon exocytosis in our model (not shown) because these 
channels have an insignificant influence on the AP upstroke (see 
above in Sec. 3.3.1). 

However, the L-type Ca 2 * channel blocker isradipine sup- 
pressed the electrical activity induced by low glucose or tol- 
butamide in isolated a-cells and glucagon secretion in human 
islets. 8 This raises the question of whether the change in activity 
of L-type Ca 2 * channels can lead to appropriate changes in AP 
firing. We found that suppression of AP firing with L-type Ca 2 * 
channel block can be obtained by changing the basal conduc- 
tivities of some channels in our model. For example, a decrease 
in conductivity of voltage-gated Na* channels and increased 
conductivity of L-type Ca 2 * channels in comparision with basal 
levels in Table 2 can lead to an increased effect of L-type Ca 2 * 
current on spike formation. In such conditions decreased con- 
ductivity of voltage-gated L-type Ca 2 * current can suppress AP 
firing (Fig. 9). Based on this evaluation we can suggest that 
human a-cells might have increased L-type Ca 2 * channel activ- 
ity in comparison with mouse a-cells. Interestingly, small spikes 
were simulated in this case after block of L-type Ca 2 * current (as 
well as in Fig. 6 at block of K Dr channels) that can also be found 
after suppression of AP firing in several experimental conditions 
(see for example 8,54 ). 

3.3.8. Non-L-type Ca 2 * channels. Non-L-type Ca 2 * channels 
include N- and P/Q-type VGCCs. These channels have been 
linked to glucagon exocytosis and may be opened at PM poten- 
tials that are above that for opening L-type Ca 2 * channels, 31 and 
we have modeled this high-voltage activation (Supp. Material, 
Eq. A28andFig. 10). 

An N-type Ca 2 * channel blocker (w-conotoxin) inhibited 
glucagon release elicited by low glucose alone in mouse islets 29 ' 39 
and in single rat a-cells. 27 P/Q-type Ca 2 * channels can serve as 
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the main component of non-L-type Ca 2 * channels in triggering 
exocytosis of glucagon granules in human a-cells. 8 This group 
found that a P/Q-type Ca 2 * channel blocker (w-agatoxin) did not 
inhibit spontaneous [Ca 2+ ] c oscillations but suppressed exocyto- 
sis in isolated human a-cells and inhibited glucagon secretion in 
isolated human islets at low glucose levels. 8 Maintenance of spon- 
taneous [Ca 2 *] c oscillations suggests that these channels do not 
have any important effect on AP firing (see below in Sec. 4.4). 
On the basis of these data we can suggest that block of P/ Q-type 
Ca 2 * channels in human a- cells leads to suppression of glucagon 
secretion but AP firing is maintained. 

We simulated the action of an non-L-type Ca 2 * channels 
blocker by reducing their maximal conductivity (Fig. 10). The 
results of simulation correspond to the suggestion that the non- 
L-type Ca 2 * channel blocker can have only weak effects on elec- 
trical activity; however it can still significantly suppress glucagon 
secretion. This is because non-L-type Ca 2 * channels play a small 
role in the AP upstroke but they are strongly responsible for glu- 
cagon release (see above in Sees. 2.4 and 3.3.1). 

3.3.9. Ct channels and y-aminobutyric acid ( GABA) action. 
GABA receptors have been identified in a-cells and glucagon 
secretion can be inhibited by GABA release from neighboring 
(3-cells in rat. 3,55 ' 56 The existence of CI" currents activated by 
GABA in a-cells was described for cells isolated from guinea 
pig and rat islets, and application of GABA led to hyperpolar- 
ization and termination of AP firing. 28,55 Interestingly, activa- 
tion of GABA A receptors induces an inward CI" current that 
inhibits neuronal firing in the CNS as consequence of PM 
hyperpolarization. 57 We have modeled the GABA-activated 
inward CI" current (Fig. 1 and Supp. Material, Eqn. A39). 
Simulation of an activation of this current (as an increase in 
its maximal conductance) led to PM hyperpolarization that 
abolished AP firing and decreased glucagon secretion. For 
example, the simulated result (at g mCIG = 40 pS and above) was 
similar to Figure 7 (not shown). This mechanism is similar 
to a PM hyperpolarization and a glucagon secretion decrease 
by an activation of K ATp channels (see for example left part of 
Fig. 5). 

3.4. Ca 2 * dynamics and glucagon secretion. Increased glu- 
cagon secretion at low glucose levels is often associated with ele- 
vated [Ca 2+ ] c in a-cells. 40,58,59 However, a glucose-induced drop in 
[Ca 2+ ] c in a-cells likely contributes only partly to the strong glu- 
cose-induced inhibition of glucagon secretion in islets. 41 [Ca 2+ ] c 
elevation was also found in a-cells of isolated mouse islets simul- 
taneously with decreased glucagon secretion following exposure 
to increased glucose. 23 

We have analyzed these contradictory experimental evidences 
using our model. Cytoplasmic Ca 2+ concentration is a result of 
the interaction of the characteristic set of Ca 2 * channels, Ca 2+ 
pumps and electrical activity. Ca 2+ enters cells through all types 
of open Ca 2+ channels. On the other hand, glucagon release 
seems to be regulated mainly by specific high-voltage-activate 
Ca 2 * channels (non-L-type, as discussed above in Sees. 2.4 and 
3.3.8). This opens up the possibility for separate dynamics of 
[Ca 2 *] c and glucagon secretion, in contrast to (3-cells (see ref. 31). 
Our simulations illustrate these cases. For example, [Ca 2 *] c and 



0 r 



I" 25 
> -50 

-75 

B 1.5 

S? 1 
O 0.5 
0 




■ ■ I ■ ■ I— 



—I I I ■ ■ ■ I 1— 



4 8 12 16 
time (sec) 



Figure 9. Simulated glucose induced spike pattern and glucagon secre- 
tion in response to the L-type Ca 2+ channel blocker at increased L-type 
Ca 2+ channel conductance (to compare with basal level) (Illustration 
to Sec. 3.3.7). (A) PM voltage and (B) relative glucagon secretion rate 
transients. The maximal conductance (g mCaL , Eqn. A22) was decreased 
from 0.4 nS to 0.05 nS at arrow. g mNa (Eqn. A14) was decreased from 32 
nS (basal level, Table 2) to 22 nS, g „„„ (Eqn. A1) was decreased from 

^mKATP ' 

basal level (30 nS, Table 2) to 25 nS and P mCap (Eqn. A38) was increased 
from basal level (300 fA, Table 2) to 500 fA to compensate for an influ- 
ence of increased L-type Ca 2+ channel current. These changes lead to 
a possibility that to AP firing may be blocked in response to the L-type 
Ca 2+ channel blocker. Other coefficients and initial conditions were as in 
Figure 3. 




Figure 10. Simulation glucose-induced spikes activity in response to 
non-L-type Ca 2+ channels blockers. (A) Action potential (V p ) and (B) 
relative glucagon secretion rate (Gs) transients. Block of non-L-type Ca 2+ 
channels was simulated by a decrease of their maximal conductance 
at low glucose level (ATP/ADP = 2). At arrow the g mCaNL (Eqn. A28) was 
decreased from 0.3 nS (basal level, Table 2) to 0.0015 nS. This leads to a 
decreased glucagon secretion at constant glucose-induced spike activ- 
ity. Other coefficients were as in Figure 3. 
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Figure 11. Simulated simultaneous [Ca 2+ ] c increase with decreased glu- 
cagon secretion due to changes in basal conductance of Ca 2+ channels 
(illustration to Sec. 3.4). (A) Action potential (V p ), (B) relative glucagon 
secretion rate and (C) [Ca 2+ ] c transients. Non-L-type Ca 2+ channels con- 
ductance (g mCaNL , Eqn. A28) was decreased from basal level (0.3 nS, Table 
2) to 0.15 nS. L-type Ca 2+ channels conductance (g mCaL , Eqn. A22) was 
increased from basal level (0.1 nS Table 1) to 0.3 nS. Glucose-induced 
changes were simulated as in Figure 3 by a step increase of the [ATP]/ 
[ADP] ratio from low (ATP/ADP = 2) to the intermediate glucose level 
(ATP/ADP = 10) at the arrow. All other parameters were taken from the 
basic set of parameters (Tables 1 and 2). 




Figure 12. Simulated glucose induced spike pattern and glucagon 
secretion at increased conductance of K ATp channels (Illustration to Sec. 
3.5). (A) Action potential (V ) and (B) relative glucagon secretion (Gs) 
rate transients. Maximal K flTp channels conductance (g mKATP in Eq. A1) was 
increased from 30 nS (basal level) to 75 nS. Glucose-changes were simu- 
lated (as in Fig. 3) by a step increase of the [ATP]/[ADP] ratio at the arrow 
from a low (ATP/ADP = 2) to an intermediate glucose level (ATP/ADP = 
10, that corresponds to 4.7 mM glucose in Fig. 2D). All other parameters 
were taken from the basic set of parameters (Tables 1 and 2). AP firing 
was induced only at high glucose levels. 



glucagon secretion have changed in one direction in the simula- 
tion of increased glucose using a chosen basal set of parameters 
(Fig. 3). However, the simulation in Figure 11 shows that PM 
depolarization following increased glucose can lead simultane- 
ously to increased [Ca 2+ ] c and decreased glucagon secretion when 
the conductivity for L-type channel was increased and the con- 
ductivity a non L-type Ca 2t channel was decreased (in compari- 
son with the basal set in Table 2). 

This result shows that when Ca 2+ enters cells mainly through 
Ca + channels which also promote glucagon release (non-L-type) 
the [Ca 2+ ] c and glucagon secretion should change in the same 
direction with increased glucose (as in Fig. 3), because the frac- 
tion of Ca 2t that penetrates into cell through L-type Ca* chan- 
nels during spikes is small. However, glucagon secretion can 
be decreased as a consequence of suppression of AP firing and 
closing non-L-type Ca + channels, but [Ca 2+ ] c can increase fol- 
lowing an increase of Ca 2+ influx through L-type Ca 2t channels 
following further PM depolarization, if the fractional Ca 2+ 
entering via L-type Ca + channels is increased. These exam- 
ples show that Ca 2+ dynamics may depend on the relationship 
between Ca + channels that have different voltage-dependent 
characteristics. 

However, [Ca 2t ] c and glucagon secretion should always 
decrease simultaneously with PM hyperpolarization that blocks 
all high voltage-gated Ca 2t channels (simulations not shown). 
This conclusion is consistent with the experimental data that 
opening of K AT[) channels with diazoxide that leads to PM repo- 
larization lowered [Ca 2t ] c in isolated mouse a-cells. 24 ' 41 

3.5. Experiments with activation of AP firing and glucagon 
secretion at high glucose levels. Isolated, purified a-cells are 
removed from any potential paracrine regulation from neigh- 
boring (3- and 8-cells and only intrinsic glucose regulation can 
take place. Two different patterns were found in the investiga- 
tions of isolated rodent a-cells. Isolated rat pancreatic a-cells 
were electrically active in the absence of glucose but glucose 
(20 mM) inhibited electrical activity. 40 Glucose reduced also AP 
firing leading to PM depolarization and an inhibition of gluca- 
gon secretion in isolated normal mouse a-cells. 24 Obviously, this 
pattern corresponds well to the metabolic mechanism of regula- 
tion of glucagon secretion (as discussed above in Sec. 3.2 and 
Fig. 3). 

Other investigators found that glucose stimulated rather than 
inhibited glucagon release in isolated rat a-cells. 27 ' 54,S0,61 Recent 
measurements of glucagon secretion from pure populations of 
flow-sorted mouse a-cells also have shown that contrary to its 
effect on intact islets, glucose stimulated glucagon secretion. 23 
Indeed, the PM was repolarized and AP firing was suppressed at low 
glucose but AP firing appeared and glucagon secretion increased 
with increased glucose in isolated rat a-cells. 54 Depolarization 
from a hyperpolarized level and an appearance of AP firing 
were also observed in response to glucose in dispersed mouse 
a-cells 11 ' 44 

We were able to simulate this pattern using our model. For 
example, we found that AP firing is absent at low glucose lev- 
els but can be stimulated with increased glucose (Fig. 12) if 
the maximal conductance of K ATp channels was increased in 
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Figure 13. Schematic diagram of metabolic and paracrine mechanisms at electrophisiologycal regulation of glucagon secretion. (A) Metabolic mecha- 
nism of regulation. This mechanism suggests that glucose metabolism leads to K ATp channel inhibition at high glucose, which depolarizes the a-cell, 
inhibiting action potential firing and glucagon release (as in Fig. 3). Glucose-changes were simulated (as in Fig. 3) by a step increase of the [ATP]/[ADP] 
ratio at the arrow. V p is PM potential, Gs is the relative rate of glucagon secretion. (B) Metabolic mechanism with additional PM hyperpolarization. 
Additional PM hyperpolarization (for example by an increase of K flTp channels conductivity, see Sec. 3.5 and Fig. 12) leads to suppression of AP firing 
and glucagon release at low glucose level. However, AP firing activates with increase of glucose as a consequence of K flTp channel inhibition. Glucose- 
changes were simulated (as in Fig. 3) by a step increase of the [ATP]/[ADP] ratio at the arrow. (C) Paracrine interaction in regulation of islet glucagon re- 
lease. At low glucose level, where paracrine influence is insignificant, intrinsic a-cell mechanisms of channel interaction lead to AP firing and increased 
glucagon release. Increased glucose level causes insulin, zinc and GABA release from (3-cells and somatostatin from 8-cells. This leads to membrane 
hyperpolarization, block of AP firing and inhibition of glucagon release (Sec. 4.2). Paracrine action leading to PM repolarisation was simulated as in 
Figure 7 at the arrow. (D) Glucagon secretion with suppressed both metabolic and paracrine regulatory mechanisms. Intrinsic a-cell mechanisms of 
channel interaction lead to AP firing and increased glucagon release independently on glucose concentration if both metabolic and paracrine regula- 
tory mechanisms are suppressed or can compensate for each other. Glucose-change is represented by arrow. 



comparison with the basal level (that we used for simulation in 
Fig. 3), even though ATP/ADP ratio changes were the same as 
describing a simulation of a-cell glucose sensing in Figure 3. 
This can be attributed to the additional hyperpolarization of 
PM with increased / KATP leading to suppression of AP firing and 
glucagon secretion at low glucose. [We have considered above 
the similar suppression of AP firing and glucagon release at the 
increase of maximal conductance of K ATp channels at action of 
activators (left part of Fig. 5)]. It is noteworthy that any influ- 
ence on ion currents leading to additional PM hyperpolariza- 
tion at low glucose levels (other than voltage-gated currents) 
can lead to a similar effect. For example, the corresponding 
somatostatin activation of the G protein-gated K* channel or 
GABA activation of CI" current (that were suggested as zero for 
basal set of coefficients) can lead to the pattern in our model 
that is similar to Figure 12 (not shown). 

These data and our previous analysis show that increased 
glucagon secretion with increased glucose from a low to high 
level can take place if AP firing was suppressed at low glu- 
cose but can be induced by PM depolarization with increased 
glucose for example as a consequence of closing K ATp chan- 
nels (see Fig. 12 and Sec. 3.5). We suggest that the metabolic 
mechanism of regulation of glucagon secretion can act in both 
cases. This leads to PM depolarization with an increase of glu- 
cose leading to different patterns that can depend on initial 
PM repolarization with an associated increase in threshold 
potential in the second case that blocks AP firing at low glu- 
cose. Such differences in PM polarization in isolated a-cells 
are difficult to explain, but may be due in part to different 



methods of a-cells isolation leading to different conductivities 
of channels. However, these differences require further 
investigations. 

4. Computational Systems Analysis of the 
Mechanisms of Glucagon Secretion in Pancreatic 
a-cells 

We analyzed the mechanism of metabolic regulation of glucagon 
secretion in Sec. 3. In this section we will analyze other mech- 
anisms for regulation of AP firing and glucagon secretion that 
were described in literature. Schematic diagrams for these dis- 
cussed mechanisms are presented in Figure 13. 

4.1. Regulation of glucagon secretion through an effect 
on electrical activity. Numerous experimental studies and our 
simulations show that glucagon secretion can be regulated by 
electrical activity in a-cells if the opening of specific VGCCs is 
required for glucagon secretion (see above in Sec. 3). The ATP- 
sensitive K + channel is considered to be the primary channel that 
can regulate PM potential. Experiments in which increasing 
concentrations of tolbutamide and diazoxide were used to titrate 
K ATp -channels activity in a-cells in intact mouse, rat and human 
islets, indicate that an increase of glucagon secretion only occurs 
within a narrow window of K ATp -channel activity where AP fir- 
ing occurs. 3 ' 8 ' 30 ' 41 We were able to simulate these results by chang- 
ing K ATp -channel activity (Fig. 5). However, our simulation also 
shows that PM depolarization can increase glucagon release even 
without activation of AP firing simply by opening additional spe- 
cific VGCCs (right part of Fig. 5). 
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In general, these experiments and our simulations suggested 
that glucagon release may depend on membrane potential in a 
complex way that can lead to opposite results for the dependence 
on an initial PM potential. For example, PM depolarization can 
lead to AP firing and increased glucagon secretion if a-cells were 
hyperpolarized previously below threshold for AP firing and a 
depolarization activates this AP firing, or if there is an increased 
PM potential from medium to high level even though AP firing 
is absent (see Fig. 5). On the other hand, a depolarization can 
block AP firing and decrease glucagon secretion if AP firing was 
pre-existent (Figs. 3 and 13A). There are many different mech- 
anisms that can influence electrical activity and consequently 
glucagon secretion and we have attempted to analyze several of 
them. 

4.2. Electrophysiological basis of paracrine regulation of 
glucagon secretion. Rodents form a characteristic a-cell mantle 
surrounding the core of (3-cells, however, a- and (3-cells show 
a more mixed distribution and the a-cells are frequently juxta- 
posed to the (3-cells within human islets. 62,63 Products released by 
(3-cells can enter the immediate proximity of the a-cells. Thus, 
the a-cells can be exposed to high interstitial concentration of 
(3- and 8-cells secretory products whenever they are secreted. 
This interplay between neighboring islet cells is a kind of "para- 
crine interaction" that may control glucagon release in a-cells. In 
this regard, in addition to a glucose, several paracrine inhibitors 
of glucagon secretion from non a-cells (i.e., insulin, zinc ions 
or GABA from (3-cells or somatostatin from 8-cells) have been 
proposed (for a review see refs. 1,3 ' 4,64 ). 

Modulation of ion channel activity would permit paracrine 
signals to have a rapid and precise effect on glucagon secre- 
tion. 65 We analyzed this proposal using our computational sys- 
tems analysis approach. For example, our analysis showed that 
somatostatin and GABA can strongly inhibit glucagon secre- 
tion in pancreatic a-cells activating the corresponding channels 
that cause hyperpolarization of the PM, suppression of electri- 
cal activity and decreased glucagon secretion (see above in Sees. 
3.3.4 and 3.3.9). 

Insulin also has been proposed as an intraislet paracrine factor 
from neighboring (3-cells that can suppress glucagon secretion 
in rodent and human islets. 1 ' 3 ' 4,8 ' 64,66 " 68 Consistent with this role, 
the insulin receptor and the insulin-receptor downstream signal- 
ing molecules are highly expressed in pancreatic a-cells (see for 
Rev. 69 ). Recent investigations have provided some insight into the 
mechanism of insulin action. Insulin transiently inhibited elec- 
trical activity and glucagon secretion in isolated rat a-cells, most 
probably through a signaling pathway that triggered membrane 
hyperpolarization. 54 However, insulin may not modulate the spe- 
cific channels as for GABA or somatostatin. Insulin can activate 
K ATp -channels in mouse a-cells, an effect that may be mediated 
via the phosphatidylinositol 3-kinase signaling pathway 70 lead- 
ing to hyperpolarization. Interestingly, insulin can also act by 
modulating Akt, a critical downstream effector of the phos- 
phatidylinositol 3-kinase, leading to a translocation of GABA A 
receptors to the cellular membrane to allow its ligand, GABA, 
also leading to membrane hyperpolarization. 71 We have ana- 
lyzed the mechanisms of inhibition of glucagon secretion by PM 



hyperpolarisation by activation of K ATp - and GABA-dependent 
CI" channels (see above in left part of Fig. 5 and Sees. 3.3.2 and 
3.3.9). 

Zinc is present in and co-crystallizes with insulin in (3-cell 
granules and it is secreted from (3-cells on exposure to high glu- 
cose, and thus could have an inhibitory effect on rat and mouse 
glucagon secretion (for a review, see refs. 54, 64 and 72). Zn 2+ 
may reversibly activate K ATp -channels concordantly reducing 
electrical activity and glucagon secretion in isolated rat a-cells. 54 
Similarly, a decrease in Zn 2t may lead to closing K ATp channels in 
isolated mouse islets. 72 The mechanism of inhibition of glucagon 
secretion as a consequence of an activation of K ATp -channels and 
PM repolarization was analyzed above (see Sec. 3.3.2 and left 
part of Fig. 5). 

Based on the preceding analysis we can explain the manner in 
which paracrine mechanisms are able to inhibit glucagon secre- 
tion with increased glucose. Negative PM potential can exist in 
a-cells at low glucose level within a range that allows AP firing 
leading to glucagon release. Increased glucose leads to increased 
pancreatic insulin, GABA and zinc secretion from (3-cells and 
somatostatin from 8-cells. These paracrine factors should lead 
to a-cell repolarization as a consequence of activation of corre- 
sponding channels. Decreased glucagon secretion occurs when 
a-cells are hyperpolarized lower than a threshold level for AP 
firing (Fig. 13C). 

4.3. Hormonal regulation. Several hormonal regulators 
were also investigated. Leptin receptor isoforms are expressed in 
mouse and human a-cells. 73 Application of leptin hyperpolarized 
the membrane potential in glucagon-secreting alphaTCl— 9 cells 
and suppressed the electrical activity induced by low glucose. 
These effects were accompanied by decreased glucagon release 
from a TCI— 9 cells and mouse islets. 73 According to these data 
and based on our analysis decreased glucagon secretion after the 
application of leptin can be explained by PM repolarization and 
suppress of AP firing as was suggested for paracrine mechanisms 
(see Sec. 4.2 and Fig. 13C). The action of leptin leading to PM 
hyperpolarization may be mediated by opening K ATp channels in 
the mouse pancreatic (3-cells. 74,75 A similar process of K ATp chan- 
nels activation by leptin may also occur in the pancreatic a-cells 
(see also 73 ). 

The incretin hormone GLP-1 is a powerful suppressor of 
glucagon secretion. The mechanism of suppression of glucagon 
secretion seems to include direct action of GLP-1 on a-cells 
and an indirect one through activation of (3- and 8-cells. 39,76,77 
Interestingly, GLP-1 can inhibit glucagon secretion without 
much of an effect on a-cell electrical activity in isolated mouse 
islets, and it was suggested that GLP-1 can selectively inhibit 
non-L-type Ca 2+ channels in mouse which are responsible for 
glucagon release. 39 According to our simulation an inhibition of 
non-L-type Ca 2+ channel can indeed lead to inhibition of gluca- 
gon secretion without suppression of AP firing (see above, Fig. 10 
and Sec. 3.3.8). 

GLP-1 receptor mRNA was not detected in isolated rat 
a-cells nor did GLP-1 stimulate glucagon release in these cells. 54 
It was proposed that the mechanism of action of GLP-1 might 
involve inhibition of the rat a-cells by somatostatin secreted from 
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neighboring somatostatin-producing 8-cells that can be activated 
by GLP-1. 77 We considered this mechanism of somatostatin 
action above (see Fig. 7). GLP-1 can also activate insulin secre- 
tion in (3-cells, 78,79 and insulin can inhibit glucagon secretion 
(see above in Sec. 4.2). However, the molecular details of GLP-1 
action on a-cells remain elusive. For example, GLP-1 might also 
directly effects regulation of glucagon granule exocytosis (as it 
occurs in (3-cells, see refs. 16 and 79). At present these mecha- 
nisms have not been investigated in a-cells and we were not able 
to analyze them using our computational approach. 

4.4. [Ca 2+ ] c handling in od-cells. A rise in the cytoplasmic 
Ca 2+ is the important physiological trigger for fusion of the 
insulin-containing granules in (3-cells (for a review see refs. 
16 and 80). However, global [Ca 2+ ] c may not directly deter- 
mine glucagon release in a-cells even though a calcium influx 
through VGCCs located in the plasma membrane can play a 
decisive role in the fusion process for exocytosis of granules 
both in a- and |3-cells (as discussed above in Sees. 2.4 and 3.4). 
We found that the difference in the voltage dependence of cor- 
responding Ca 2t channels can explain contradictory data about 
cytoplasmic Ca 2+ and glucagon secretion dynamics (see above 
in Sec. 3.4). 

Isolated a-cells in intact mouse and human islets exposed to 
low glucose (or in the absence of glucose) exhibited spontane- 
ous [Ca 2+ ] c oscillations with period up to several minutes (for 
a review see refs. 8 ' 9,23,39 ' 41,59 )- Spontaneous [Ca 2+ ] c oscillations in 
isolated mouse a-cells result from synchronous spiking electri- 
cal activity. 41 The mechanisms that generate and terminate AP 
spikes are responsible for the periodical AP firing (bursting) in 
(3-cells (see ref. 18). The periodical AP firing in (3-cells provides 
an effective mechanism for rapid changes in extracellular Ca 2+ 
influx through VDCCs that leads to cytosolic [Ca 2+ ] c oscilla- 
tions. 14,81,82 For these reasons, we also suggest that periodic AP 
firing (bursting) can also result in [Ca 2t ] c oscillations in a-cells. 
Then, an occurrence of [Ca 2+ ] c oscillations reflects the existence 
of periodic AP firing. However, our analysis has shown that block 
of AP firing can occur both during PM repolarization and depo- 
larization in a-cells (see discussion above in Sec. 3.3), and it may 
be difficult to recognize the direction of PM potential change 
(i.e., repolarization or depolarization) if only the block of [Ca 2t ] 
oscillations is detected. 

c 

Subsequent generation and termination of AP firing may be 
determined by small cyclic changes in some voltage-gated and/or 
voltage independent gating variable currents in (3-cells. 1418,83 We 
have recently performed a detailed analysis of possible mecha- 
nisms of oscillatory processes in (3-cells. 18 In light of this analysis, 
the mechanisms of [Ca 2+ ] c oscillations and periodical hormonal 
secretion in the isolated a-cells compared with (3-cells could be 
more closely examined in the future. 

4.5. Regulation of exocytosis of glucagon-containing gran- 
ules. Glucagon is secreted by exocytosis of large dense-core 
vesicles. Regulation of glucagon-containing granule secretion 
in a-cells is both complex and multifactorial and not limited to 
the effect of electrical activity. 1,5,34,84 Regulation of insulin— con- 
taining granule secretion in (3-cells is also very complex, but the 
hierarchy and interplay between the different levels of control 



are understood well enough to form a broad outline (for a review 
see refs. 79 and 80). Mathematical modeling of the regulating 
pathway and its influence on insulin granule secretion in (3-cells 
was discussed. 16,85 Similarly, several proposals have been made for 
the influence of paracrine and endocrine factors on second mes- 
sengers pathways (such as cAMP) and on glucagon-containing 
granule secretion, 1,86,87 and a kinetic model for glucagon secretion 
that describes Ca 2+ triggered exocytosis in a-cells was recently 
developed. 59 However, in contrast to (3-cells, our knowledge of 
the molecular regulation of exocytosis of glucagon-containing 
granules remains fragmentary. For this reason, we were not able 
to apply directly our computational approach to these mecha- 
nisms of regulation of glucagon release. 

4.6. a-cells modeling approach. There are a large number 
of studies on computational modeling of (3-cell function (for 
review see refs. 14, 15, 18 and 88). However, we have found only 
two computational theoretical models in the literature that con- 
sider intracellular regulatory electrophysiological mechanisms in 
a-cells. In one, a Hodgkin-Huxley-type ionic model for action 
potentials in a-cells was developed. 89 Simulation of AP firing was 
obtained in this model and it was shown that block of K ATp chan- 
nels can lead to a suppression of AP firing. However, no attempts 
were made for a simulation or analysis of metabolite changes, 
paracrine regulation, [Ca 2+ ] c handling and glucagon secretion. 
Ion channels involved in the electrical activity in mouse a-cells 
were modeled using Monte Carlo algorithms to fit steady-state 
channel currents. 90 However, simulation of AP firing or glucagon 
secretion was not considered. 

4.7. Mechanisms of regulation of glucagon secretion in 
non-human mammals. There is a significant difference in the 
mechanisms of metabolic and paracrine inhibition of glucagon 
secretion following increased glucose. According to our analysis 
both of these mechanisms of inhibition are based on a suppres- 
sion of AP firing. However, PM depolarization should be initiated 
as an effect of the metabolic mechanism, but PM repolariza- 
tion could play a decisive role for the paracrine and endocrine 
factors. This raises the question as to what mechanisms act in 
islets? 

In this light, knowledge of both electrical activity and gluca- 
gon secretion is required to establish the mechanism for stimu- 
lus-secretion coupling in a-cell. However, electrical activity of 
a-cells was not measured in many studies performed on isolated 
islets so that analysis is restricted to only a few papers. 

4.7.1. Isolated rodent islets. Glucose inhibited glucagon secre- 
tion in isolated mouse and rat islets. 3,27,39,76 However, we were able 
to find only one article where AP firing was measured directly in 
a-cells inside rodent islets. PM depolarization and AP firing were 
found in a-cells in mouse and rat isolated islets after a change 
of glucose concentration from 11.1 to 2.8 mM and hyperpo- 
larization and suppression of AP firing after return to 11.1 mm 
glucose. 12 Interestingly, other investigators 91 followed changes in 
membrane potential of individual cell types in whole mouse islets 
using a fluorescent membrane potential sensitive probe. They 
demonstrated that an increase in glucose concentration from 3 to 
10 mM is paralleled by repolarization of a-cells and depolariza- 
tion of (3-cells. Our consideration suggests that the registration 
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of PM repolarization with high glucose levels corresponds to the 
action of the paracrine mechanism of regulation (see Fig. 13C 
and Sec. 4.2) (in compliance with suggestion in ref. 12). 

Paracrine concept is difficult to reconcile for every condition, 
since several experiments have shown that the inhibition of gluca- 
gon secretion in isolated mouse islets already occurred at low glu- 
cose concentrations that have no effect on insulin or somatostatin 
secretion. 3 In this case the metabolic mechanism of regulation of 
glucagon secretion can explain these observations. However, AP 
firing was not measured directly in a-cells inside islets and this 
makes an explanation of these results problematic. 

4.7.2. Diabetic animal models. Alloxan and streptozotozin 
treated rodent, dog and pig islets have been studied as islet injury 
or /|3 cell depletion models. When glucose was infused without 
insulin glucagon did not decrease but rose up to 50% in a dia- 
betic dog. However, when large amounts of insulin was infused 
together with glucose, a prompt decline in glucagon (4-fold) 
was observed. 92 Postprandial insulin-driven suppression of glu- 
cagon secretion was also lost in the pigs treated with alloxan. 93 
Similarly, glucagon secretion did not increase during hypoglyce- 
mia in streptozotocin treated rats, but it could be increased if the 
intrapancreatic insulin was switched off during hypoglycaemia. 94 
A decisive role for endogenous insulin for the glucose-induced 
suppression of mouse glucagon secretion was found for strepto- 
zotozin treated mice. 68 Obviously, a metabolic mechanism could 
not explain the suppression of glucagon secretion in these cases 
but a paracrine mechanism may do so. 

Importantly, our analysis of experimental observations 
obtained in vivo leads to the conclusion that paracrine mecha- 
nisms acting through PM repolarization may play a decisive role 
in suppression of glucagon secretion following increased glucose 
in non-human mammals. 

Other proposals have also been made concerning mechanisms 
of regulation of glucagon release (such as direct effects on gran- 
ule exocytosis, see above in Sec. 4.5) besides mechanisms act- 
ing through changes of PM potential. However, these proposals 
require further investigations and at present they are beyond the 
scope of our computational systems analysis. 

4.8. Comparison of the mechanisms of regulation in a-and 
(3-cells. It may seem puzzling that a-cells, having a similar met- 
abolic pathways and ion channel set as the (3-cells (at least in 
human, see refs. 8 and 95) can be regulated oppositely to (3-cells. 
It was proposed that differences in ion channel expression con- 
tribute to the opposing regulation of insulin and glucagon secre- 
tion following changing levels of glucose. 2 Our computational 
analysis indicates how the variations in the conductivities of the 
channels can lead to a switch from characteristic a-cell AP firing 
at low glucose levels to AP firing only at increased glucose. We 
have shown in silico that an activation of the relatively strong 
hyperpolarizating currents in a-cells can perform such a shift 
(see above in Sec. 3.5). For example, the increase in K ATp channel 
conductance is enough to transform AP firing from a- to (3-cell 
type behavior (Fig. 12). Interestingly, a relatively higher density 
of K ATp channels in (3- cells than a-cells was reported. 11 ' 44 ' 96 ' 97 
According to our analysis this can be sufficient to explain the 
differences in glucose dependent electrophysiological properties 



a- and (3-cells. However, any influence on ion currents resulting 
in additional PM hyperpolarization at low glucose levels may lead 
to a similar effect (see Sec. 3.5.). For this reason, it is also possible 
that (3-cells have decreased activity of Na* background current or 
some other depolarizing factors in comparison with a-cells that 
also can lead to a corresponding shift in AP firing in region of 
increased glucose. 

Another important difference is the possibility that over-all 
cytoplasmic calcium concentration may not play a decisive role 
in an activation of glucagon release in a-cells (as compared with 
microdomain sequestered calcium entry) whereas insulin secre- 
tion depends significantly on global cytoplasmic calcium (see 
discussion above in Sec. 3.4). 

PM depolarization as well as a hyperpolarization can lead to 
block of AP firing and associated decrease of glucagon secretion 
in a-cells (see above in Sees. 3.3.2 and 4.1). This is a significant 
distinction from (3-cells where only hyperpolarization leads to a 
decrease of insulin secretion in physiological conditions. 

The molecular characterization of a metabolic glucose sensor 
is most worked out in pancreatic (3-cells, with important roles 
for glucokinase and mitochondrial oxidative fluxes in the regula- 
tion of ATP-sensitive K + channels. We suggest that the metabolic 
mechanism of PM depolarization by K ATp channel inhibition 
with increased ATP/ADP ratio is significantly suppressed in 
a-cells (see above in Sec. 3.1). Specific components of this meta- 
bolic mechanism which are well expressed in (3-cells (such as a 
low activity of lactate dehydrogenase and uncoupling protein 2 
and increased activity of pyruvate dehydrogenase) are compro- 
mised in a-cells (see above in Sees. 2.1 and 3.1). In this case the 
metabolic mechanism can serve to maintain a sufficient ATP/ 
ADP ratio at low glucose levels that leads to regulation of the 
PM potential within narrow limits where AP firing can occur. 
However, the increase in ATP/ADP ratio following a further 
increase in glucose is supressed in a-cells (see discussion above 
in Sec. 3.1). An important component of metabolic glucose sens- 
ing such as K ATp channels may play an important role in regula- 
tion of a-cell function through a change in PM potential. For 
example, a repolarization of PM through a paracrine mechanism 
action may occur through opening of K ATp channels (see Sec. 
4.2). 

We suggest that such a-cell behavior is not surprising, if an 
aim of this system is a suppression of glucagon secretion follow- 
ing increased glucose. Indeed PM depolarization is the mecha- 
nism leading to insulin secretion in (3— cells. For this reason 
we suggest that PM repolarization should be responsible for 
decreased glucagon secretion that according to our analysis may 
be a consequence of paracrine mechanisms. Block of a metabolic 
mechanism can open the way for this kind of mechanism for PM 
repolarization (see above in Sec. 4.2). 

We hypothesize that the suppression of a glucose sensitive 
mechanism in the a-cell can also lead to different fates of a- and 
(3-cells during the developing of type 2 diabetes mellitus. The 
existence of specific mechanisms for (3-cell glucose sensitivity 
itself, and particularly the increase of mitochondrial membrane 
potential with increased glucose, may underlie the increased sen- 
sitivity of these cells to injury (see for details, ref. 25). According 
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to this hypothesis the absence of a significant increase of mito- 
chondrial potential with glucose, as simulated in our model (Fig. 
2B), may prevent a-cells from metabolic injury, keeping them 
intact when B-cells can be damaged. 

5. Human a-cell Function and Diabetes 

The growing intense focus on the role of the a-cell in diabetes 
(for a review see refs. 6, 19 and 20) has increased the importance 
of understanding how electrical activity induced by metabolism 
of sugar and other agonists is established and how it is altered in 
various pathological states. 

5.1. Human isolated a-cells and islets. Human islets contain 
a substantially higher proportion of a-cells than rodent islets, 
normal a-cell mass in the whole human pancreas amounting to 
-40% of B-cell mass. 58 Comprehensive electrophysiological char- 
acterization of the ion channels was made recently for isolated 
human a-cells. 8 The same set of channels exist in mouse and 
human a-cells (see refs. 8 and 31). 

Regenerative electrical activity was observed in the major- 
ity of isolated human a-cells analyzed at low glucose levels. 
Tolbutamide (100-400 |xmol/i) depolarized the cells, and was 
associated with a reduction of AP peak voltage. 8 These results 
can be explained by some role of K ATp channels as a factor that 
regulates PM potential similarly to isolated mouse a-cells (see 
Sec. 3.3.2). However, electrical activity of a-cells inside human 
islets was not measured. 

Intracellular Ca 2+ oscillations in a-cells inside isolated human 
islets were registered at low glucose levels and these oscillations 
were suppressed with increased glucose." This pattern can corre- 
spond to existence of AP firing at low glucose and its suppression 
with increased glucose level (see discussion above in Sec. 4.3). 

Increased glucose concentration inhibited significantly gluca- 
gon secretion at 6 mM in comparison with 1 mM, and insu- 
lin inhibited glucagon secretion at low glucose level (ImM) in 
isolated human islets. 8 Insulin secretion increased also sharply 
in isolated human islets (other than in isolated mouse islets at 
increased glucose level from 1 to 5 mM. 3,8 Inhibition of GABA A 
receptor increased glucagon secretion significantly at high glu- 
cose levels, however this inhibition did not effect insulin secre- 
tion in isolated human islets. 56 These data can be interpreted 
as evidence of suppression of glucagon secretion by insulin and 
GABA at elevated glucose. B-cell-mediated paracrine signaling 
may predominate over intrinsic signals in the regulation of gluca- 
gon secretion in human and there are numerous indirect data and 
speculations that intrinsic insulin can be major paracrine factor 
for a-cells in human islets. 4,6 ' 64 

Distinct paracrine factors may play a different role in rodent 
and human islets. For example, it was found that Zn 2+ stimu- 
lated rather than inhibited glucagon secretion in isolated human 
islets. 8 Zn 2t did not play a significant role in the regulation of glu- 
cagon release in type 1 diabetic humans where only a few B-cells 
exist. 100 

Based on these limited date and our previous consideration we 
can suggest that a-cells in normal human islets (as well as in non- 
human mammals) have AP firing at low glucose levels with an 



associated increase in glucagon secretion. AP firing and glucagon 
secretion are suppressed by paracrine factors (possible by insulin) 
with increased glucose and a metabolic regulatory mechanism 
may not play a significant role. However, further studies are nec- 
essary to investigate this important question. 

5.2. Abnormal a-cells function and diabetes. Normally, the 
postprandial increase in plasma glucose leads to increased insulin 
secretion and the moderate suppression of glucagon release. On 
the other hand, the physiological defense against falling plasma 
glucose includes decreased pancreatic insulin, Zn 2+ and GABA 
secretion from the B-cell and somatostatin from 8-cells, as well 
as an increased glucagon secretion that prevent a sharp decrease 
of glucose concentration in blood. 6,20 

These defense mechanisms are compromised in type 1 and 
type 2 diabetes where B-cell function is significantly inhibited. 
Glucose sensing in glucagon release becomes defective. Diabetes 
is associated with aberrant secretion of glucagon at high plasma 
glucose levels (normally inhibited at such levels). 6,19,101,102 In 
patients with type 2 diabetes it has been shown that the increased 
glucagon concentration postprandially contributes to hypergly- 
cemia because hyperglucagonemia results in increased hepatic 
glucose production. This adds an endogenous source of glucose 
to the exogenous glucose from the meal. 103,104 However, mean 
fasting plasma glucagon concentration may not differ signifi- 
cantly between nondiabetic and type 2 (often and type 1) dia- 
betic subjects. 6,102,105 

Although the precise mechanisms of relative hypergluca- 
gonemia with increased plasma glucose level in the diabetic 
state is still obscure, this behavior can be caused by loss of the 
decrement in intraislet paracrinic factors from B-cells (insu- 
lin, GABA and/or zinc) that normally results in suppressing of 
glucagon secretion at high glucose level (as discussed above in 
Sec. 4.7.2). 

AP firing is required to maintain a significant rate of glu- 
cagon secretion (see discussion above in Sees. 3.2 and 5.1). 
Maintenance of AP firing and the associated high rate of gluca- 
gon release at elevated glucose level can readily be obtained in 
our model by blocking paracrine mechanisms of glucose sens- 
ing (Fig. 13D). It is also possible that the metabolic mechanism 
leading to PM depolarization under conditions of low paracrine 
activity, that usually leads to hyperpolarization, can compensate 
for each other at high glucose levels to continue AP firing. 

Interestingly, elevated plasma glucagon levels are associated 
in type 1 diabetes with severe diabetic ketoacidosis. 102,106 While 
mechanisms underlying hyperglucagonemia in these cases are 
still obscure this may be due to the combination of low insu- 
lin concentrations and B-cell contents in the islet. For example, 
insulin can regulate gene expression in a-cells through several 
transcription factors (see ref. 69). 

In islets isolated from donors with type-2 diabetes insulin 
release was increased (however significantly less than for non- 
diabetic donors) and glucagon release was decreased (however 
this decrease was less than for nondiabetic donors) at transitions 
from low (1 mM) to high (16.7 mM) glucose level. 56 These data 
correspond to an expected decrease in paracrine inhibition glu- 
cagon release with decreased activity of B-cells. 
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However, glucagon release from islets isolated from donors 
with type-2 diabetes was reduced at low glucose and that glucose 
significantly stimulated rather than inhibited secretion in these 
islets in work. 3 Based on our foregoing analysis we can explain 
this observation. For example, at suppressed paracrine mecha- 
nisms this behavior resembles the second case for isolated rodent 
a-cells (see Fig. 13B and Sec. 3.5) and can be explained in the 
same manner, that is the initial PM potential could be lower than 
the threshold for AP firing in these islets and AP firing is then 
activated when glucose is increased leading to PM depolarization. 
However, further studies, including electrophysiological record- 
ings with secretion measurements, are necessary to determine 
the mechanisms of regulation of glucagon secretion in diabetic 
human islets. 

5.3. Potential therapeutic approaches. New therapeutic 
approaches targeting glucagon secretion may prove useful in the 
treatment of diabetes. 19,20 According to our analysis the paracrine 
mechanism acting through PM repolarization may play a decisive 
role in a suppression of glucagon secretion following increased 
glucose in human as well as in non-human mammals. For this 
reason we can suggest that agents inducing PM hyperpolarization 
with increased glucose may be useful as a therapeutic approach. 
For example, somatostatin infusion inducing PM repolarization 
(see above in Sec. 4.2) was used to suppress hyperglucagonemia 
in type 1 diabetes. However, several significant side effects of 
somatostatin and its long acting analogs preclude its chronic use 
(see refs. 20 and 107). 

Leptin can also hyperpolarize the PM in a-cells (see above 
in Sec. 4.3). For this reason, the inhibitory effect of leptin on 
the a-cells may be a plausible therapeutic strategy for suppression 
of glucagon secretion in Type 1 diabetes. 20,108 However, leptin 
could also suppress insulin secretion 74 so this strategy may have 
complex effects on both glucagon and insulin secretion if 3-cells 
maintain activity (see also refs. 20,73). 

Amylin (IAPP, or islet-ssociated polypeptide) is a (3-cell hor- 
mone that is co-secreted with insulin. Clinical studies, in which 
a commercially available amylin analog (pramlintide) was used, 
have shown that amylin can decrease circulating glucagon lev- 
els. 20,109 The mechanism of amylin analog action, as well as other 
paracrine factors, may lead to PM hyperpolarization in a-cells. 
However, further studies, including electrophysiological investi- 
gations, are necessary. 

Interestingly, K ATp channels openers (e.g., diazoxide) can 
hyperpolarize the a-cell PM, prevent AP firing and blocking 
glucagon secretion (as discussed above in Sees. 3.3.2 and 4.1). 
However, K ATp channel activation will also suppress insulin 
secretion. Application of a K ATp channel opener before or with 
meals may decrease both glucagon secretion and hepatic glu- 
cose production in Type 1 and advanced type 2 diabetes, where 
insulin secretion is insignificant. We suggest that K ATp channels 
openers can also be potential therapeutic agents in treatment of 
diabetes. However, this has not been directly addressed in clini- 
cal studies. 

Finding specific agents that lead to PM hyperpolarization only 
in a-cells rather than (3-cells may be helpful for treatment of type 
2 diabetes where intact (3-cells remain active. This can include 



specific activators of K + and CI" outward currents or blockers of 
inward Na* current in a-cells. Paracrine factors may also lead to 
a-cell repolarization as a consequence of activation of hyperpo- 
larizing ion channels (see discussion above in Sec. 4.2) and an 
activation of specific receptors can increase the associated cur- 
rents even though the amount of paracrine factors from (3-cells 
may be greatly reduced in diabetes. 

PM depolarization in a-cells can suppress AP firing and glu- 
cagon secretion (see above in Sec. 3.2). Sulphonylurea drugs (like 
tolbutamide and glibenclamide) block K ATp -channels leading to 
depolarization. However, sulphonylurea therapy has in fact been 
reported to stimulate glucagon secretion in type-1 diabetes 105,110,111 
and in contrast to inhibit glucagon secretion in both healthy 66,105 
and type-2 diabetic subjects. 112 Our model suggests that the 
effect of sulphonylureas depends on initial PM potential condi- 
tions (see Sec. 4.1) and, in principle, can even increase glucagon 
release if AP firing is initially suppressed and PM depolarization 
leads to initiation of APs as shown in Figure 13B. Increased PM 
depolarization from a medium level (in medium glucose) even 
at suppressed AP firing can also increase glucagon secretion (see 
discussion above in Sec. 3.2 and right part of Fig. 3). This may 
explain the contradictory data on the effect of sulphonylureas on 
glucagon secretion. 

The incretin hormone GLP-1 is a powerful suppressor of glu- 
cagon secretion (see above in Sec. 4.3). GLP-1 based therapies 
target both insulin and glucagon secretion, and inhibition of 
glucagon accounts for a large part of the glucose-lowering effect 
in type 2 diabetic patients. 19,77 There is reason to believe that 
indirect reciprocal (3-cell-mediated paracrine signaling (partially 
the GLP-1 mediated increase in insulin secretion) may predomi- 
nate in the suppression of glucagon secretion in type 2 diabetic 
patients. 105 Interestingly, the use of a GLP-1 analog, liraglutide, 
may help in controlling hyperglycemia with a concomitant 
reduction in insulin dose in type 1 patients 113 indicating a poten- 
tial action of GLP-1 analogs directly on a-cells. However, the 
mechanisms behind the possible direct inhibitory effect of GLP-1 
on glucagon granule exocytosis have not been elucidated and at 
present additional studies are needed to apply the computational 
approach to this effect. 

Conclusion 

Here we have reviewed intrinsic glucose-sensing, metabolic, 
paracrine and endocrine mechanisms that control glucagon 
secretion in a-cells and proposed a computational model for 
these processes. This model was developed to gain further 
insights into coupling of metabolic processes, PM potential, 
[Ca 2+ ] c handling and glucagon release in a-cells, that remains 
incompletely understood. We have shown seemingly contradic- 
tory experimental data in this field may be synthesized into a 
coherent mathematical model. In this model the metabolic fea- 
tures and major ionic currents can be fitted to recent experi- 
mental data and expression experiments. The key feature of the 
model is that it can be used to explore mechanisms underlying 
AP spikes that regulate glucagon secretion. The model repro- 
duces a number of experimental observations ranging from 
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metabolic changes and voltage-clamp current traces to Ca 2t and 
glucagon secretion dynamics. This model can be used as the 
groundwork for in silico examination of the regulation of glu- 
cagon secretion through a change in the metabolic processes, 
channels function and [Ca 2t ] c dynamics. Understanding the 
cellular and molecular regulation of glucagon release can lead to 
refinement of current therapies and to the discovery of new tar- 
gets that improve treatment of both major types of diabetes and 
also give additional insight in several types of monogenic diabe- 
tes. However, further studies supported by improved mathemat- 
ical models are necessary to understand the stimulus-secretion 
coupling in the ot-cell. 
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